Version Date: 1 February 2008 



Quantum transport calculations for quantum 
cascade laser structures 



S.-C. Lee and A. Wacker 

Institut fiir Theoretische Physik, Technische Universitdt Berlin, 10623 Berlin, 

Germany 



Abstract 

We apply a quantum transport theory based on nonequilibrium Green's functions to 
quantum cascade laser (QCL) structures, treating simultaneously the transmission 
through the injector regions and the relaxation due to scattering in the active re- 
gion. The quantum kinetic equations are solved self-consistently using self-energies 
for interface roughness and phonon scattering processes within the self-consistent 
Born approximation. In this way, we obtain the current density J, and the average 
electronic distribution f{E) at a given energy E, as a function of applied bias. As 
a test case, we apply the theory to a GaAs/Al^^Gai-ajAs QCL structure reported in 
the literature. The theoretical results reproduce well reported voltage-current (V-I) 
measurements, and also demonstrate a population inversion at a bias that agrees 
well with the range of currents and fields at which lasing is observed. 
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Quantum cascade laser (QCL) heterostructurcs[l] arc unipolar devices 
which operate on the basis of light-emitting intersubband (intraband) transi- 
tions. They are designed with a complicated sequence of semiconductor layers 
with varying widths. This sequence of layers is repeated up to, for example, 
30 periods. Part of the semiconductor layer sequence in each period operates 
as the active region within which the lasing transition occurs. The remaining 
layers in the period play the role of an injector for the active region, as well 
as a collector of carriers from the preceding period. Fig. 1 shows an example 
of the conduction band hne-up of a QCL structure taken from Ref . [2] . 

The transport of carriers through QCL structures is governed by a com- 
plex interplay between transmission through the injector region, and relax- 
ation in the active region through scattering processes and light emission ([3], 
[4]). To investigate the carrier transport in QCL structures, we apply a quan- 
tum transport theory based on nonequilibrium Green's functions [5], which 
incorporates both transmission in the injector region and scattering-induced 
relaxation processes. The quantum transport equations for the nonequilibrium 
Green's function G'^'^^(k, E), and the correlation function G^(k, E), are solved 
self-consistently with self-energies for phonon and interface roughness scatter- 
ing processes, within the self-consistent Born approximation. Prom G"^^*(k, E) 
and G^(k, i?), we can determine the current density J, and the average elec- 
tronic distribution f{E), which gives the average occupation of the levels at 
a given energy E. 

To model a QCL structure we treat the system as a periodic superlattice 
structure where each period consists of N barriers and wells with varying 
widths. The Hamiltonian for this structure may be separated into two parts, 
H = Ho + ^^scatt- The Hamiltonian Ho contains the superlattice potential 
and, if present, a static electric field S applied in the growth direction, i.e.. 
Ho = i^sL + H£. A plane- wave basis is used in the plane of the semiconductor 
layers. In the growth direction, we use a localised, Wannier state basis [6]. 
Hence, we write 



where the index n labels a period in the superlattice, and the index v labels a 
Wannier level ip'^{z) within a period, ^^^^(k) and a^(k) are creation and annihi- 
lation operators for an electron with the two-dimensional in-plane wavevector 
k, in the z/th Wannier level, in period n. The Wannier states arc not eigen- 
states of HsL, and therefore Tf represents the off-diagonal coupling between 
Wannier levels in different periods, and E" represents the diagonal elements 
of HsL in this basis. We only take into account couplings between adjacent 
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periods. The Hamiltonian Hs, due to the electric field, is written as 



He ^ Yl / ^^'^ [ - e^^oX^k) an(k) 

n ij II 

-e£:i?r( <li(k) a^(k) + <t(k) aC+i(k) )] , (2) 



where Rf'^ = J dz ■ip^{z — Id) zip'^{z). d is the length of one period. The in- 
terface roughness and phonon scattering processes expressed by ^^scatt are 
treated in the form of self-energies which are described below. We have ne- 
glected electron-electron scattering processes in this present work because of 
the increased complexity in the theory when they are included. These pro- 
cesses should, however, be treated in a more comprehensive study at a later 
stage. 

The quantum transport equations derived from the Hamiltonian H have 
the following form ([5], [6]). The retarded Green's function G]^^^^{k, E) is given 
by the Dyson equation 

E GT,a, (k, E)-Y: [{Ho)a,0 + E^«%(k, E)] (k, E) = S^,^, , (3) 
P 

where ai and a2 are general indices that include both the period and Wannier 
level indices, e.g. cci = (m, //). The correlation function G'^^^^{'k, E) is obtained 
from the Keldysh relation: 

G^a.a.(k, ^) = E G^at/3(k, E) E<^,(k, E) (k, E), (4) 
HP' 



where ^^'^^(k,^) = [G"^^*(k, The self-energies E'^^*(k,^) and S<(k,E) 
describe the scattering processes arising from -ffgcattj and are evaluated within 
the self-consistent Born approximation. Assuming that the diagonal parts of 
G(k, £■) and S(k, i?) dominate, the self-energy for interface roughness scat- 
tering has the form: 

S<'™r;(k, ^) = E (IK.S(k - kOr) G<,,„.(k', E). (5) 



V^'^^^^(k' — k) is the matrix element for interface roughness scattering. The 
equation for S'"'^*'™"^^ is obtained by the replacements: E<'''°"g'^ ^ j]rct,rough^ 
and G'^ G^^^. For optical or acoustic phonon scattering, the self-energies 
are: 



sssr(k,£^) = E iK^^s.(k,k 
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X <;[/s(^phon) + i]G': 

■yret 



ret / 

ni',ni'\ 



k',E — -Ephon) 
+ /^(-E'phon) G'5:f*„j,(k', E + -Ephon) 
k.\E — -Ephon) — 



+ 



G, 



< 



and 



k',E + E^^on) (6) 



fB{Ephon) G^j,^„j,(k', E - -Ephon) 
+ [/^(E^phon) + 1] C^^^^^Ck', E + -Bphon)|, 



(7) 



where \/p^°° represents the interaction with optical or acoustic phonons, and 

-E'phon represents the energy of the optical or acoustic phonon. fsiE) = 1 / [exp{E / ksT) - 

1] is the equilibrium phonon distribution at a temperature T. For practical pur- 

phon / 

poses, we use momentum independent scattering matrix elements Vm°i^,nu (ktyp, k'typ), 
employing a representative momentum ktyp. We find the Green's functions C*** 
and by solving the quantum transport equations, Eqs. (3) and (4), self- 
consistently with the equations for the self-energies, Eqs. (5) - (7). Periodic 
boundary conditions are used, i.e. Gm^^nAK E) = G (rn+i)iJ,,{n+i)v{k, E - le£d). 

Once C''* and have been determined, we can evaluate experimentally- 
related quantities such as the current density J = ■^coh^+'^scatt between periods 
and 1 with 



70^1 _ 
"^coh ~ 



2epo 

n 



jdE, l^nY.iHo)^^Gt,o,{KE)h 



(8) 



and 



^°rti ~] dE, j ^-§^^[Y. Sotff (k, E) G<,,o,(k, E) 



+ So<ioi(k,^)G£V(k,^)}, (9) 

where po — nie/'KTi^- S^-*^^ is given by Eqs. (5) - (7) where the sum over n is 
restricted to n = 1. J^oh is the coherent contribution to the current density 
arising from the off-diagonal couplings in i/^, i.e., in i^gL and —e£Ri^ in 
He- Jscatt represents the current contribution due to the scattering processes 
described by Hscaxx- We note here that the division between coherent and 
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scattering contributions to the current is artificial since it depends on the 
basis set used in the calculation. If, for instance, the eigenstates of Ho were 
chosen as the basis, there would be no coherent contribution to the current. 

We apply the theory to the GaAs/Ala;Gai_a;As QCL structure (Fig. 1) 
described by Sirtori et al. in Ref. [2]. To calculate the voltage-current (V- 
I) characteristic, we evaluate the current density J. Fig. 2 shows a plot of 
voltage vs. current calculated for a lattice temperature of 77 K. A comparison 
with experimental measurements (Fig. 3(a) in Ref. [2]) shows good agreement 
between the theoretical and experimental results. For instance, at a current of 
4 A, the measured voltage is around 11 V, while the theoretical curve in Fig. 2 
shows a voltage of ~ 7 V which agrees very well with the experimental results 
when a residual series resistance of ~ 1 f2 (see Ref. [2]), probably originating 
from the cladding layers, is taken into account. We comment here on the 
finding of lotti et al. [7] , of a strong influence of electron-electron scattering 
on the I-V curve. Our calculations include interface roughness scattering, but 
not electron-electron scattering. The converse is true in the calculations of lotti 
et al. where electron-electron scattering is included, but interface roughness 
scattering neglected. In both cases, I-V curves in fairly good agreement with 
experimental measurements are obtained. We therefore conclude that interface 
roughness and electron-electron scattering play similar roles in determining 
the I-V curve, and that the I-V characteristic is not sensitive to which type of 
scattering is present in the calculation. 

In Fig. 3, we show the density of states Pdos = "2 J^i^m'^I^Ou ouO^'' ^)] 
the electron distribution function f{E) = J2iy,k'^[GQ^Q^{k, -E')]/pdos- Curves at 
three different applied voltages are shown: 0, 3, and 7 V. At V, we see, as 
expected, a thermal Fermi distribution corresponding to a lattice temperature 
of 77 K. As the applied voltage increases, the occupation probability at higher 
energies increases. In addition, phonon replicas become visible in f{E). At 
7 V, the occupation probability is higher at an energy of 0.2 eV than at lower 
energies between to 0.1 eV. This demonstrates the onset of a population 
inversion. At this voltage, the current is 4.5 A, corresponding to a current 
density of 7.5 kAcm"^. This is in good agreement with the threshold current 
density of 7.2 kAcm~^ for laser emission which is reported in Ref. [2]. 

In summary, we have applied a quantum transport theory based on nonequi- 
hbrium Green's functions to QCL structures. The initial results for transport 
properties such as current densities are in good agreement with experimental 
data reported in the literature. The calculated electron energy distribution 
also demonstrates a population inversion at a current density which agrees 
well with the range of fields and currents at which lasing is observed. 
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FIGURE CAPTIONS 



Fig. 1 Example of the conduction band line-up in a QCL structure taken 
from Ref . [2] . Each period contains 16 semiconductor layers giving a period 
length d = 45.3 nm. There are 30 periods in this structure. The in-plane 
area is 6 x 10~^ cm~^. The Wannier functions used as the basis set for the 
calculations arc also shown. 

Fig. 2 Calculated V-1 characteristic for the GaAs/Ala;Gai_3;As QCL structure 
described in Ref. [2]. Sample area 6 x 10~^ cm^. 

Fig. 3 The lower row shows the distribution function of electrons for different 
applied bias. The current and applied bias are shown with each curve. The 
corresponding electric fields are 0, 22, and 53 kV/cm, and the current den- 
sities are 0, 85, and 7500 A/cm^. The upper row shows the corresponding 
density of states Pdos (scaled with po) for each bias. 
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Fig. 3 




10 



